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Abstract 

The aim of this paper is to provide conditions which ensure that the 
affinely transformed partial sums of a strictly stationary process converge 
in distribution to an infinite variance stable distribution. Conditions for 
this convergence to hold are known in the literature. However, most of 
these results are qualitative in the sense that the parameters of the limit 
distribution are expressed in terms of some limiting point process. In this 
paper we will be able to determine the parameters of the limiting stable 
distribution in terms of some tail characteristics of the underlying sta- 
tionary sequence. We will apply our results to some standard time series 
models, including the GARCH(1, 1) process and its squares, the stochastic 
volatility models and solutions to stochastic recurrence equations. 

Keywords stationary sequence stable limit distribution weak convergence mix- 
ing weak dependence characteristic function regular variation GARCH stochas- 
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1 Introduction 

Whereas there exists a vast amount of papers and books on the limit theory for 
sums Sn = Xi + ■ ■ ■ + Xn of finite variance strictly stationary sequences {Xt), 
less attention has been given to the case of sums of infinite variance stationary 
sequences. Following classical work (for example, Gnedenko and Kolmogorov 



\2i\\. Feller (1261 ) . Petrov (j46f )). wc know that an iid sequence [Xt) satisfies the 
limit relation 

a-\S^~b,,)AY^, (1.1) 

for suitable constants a„ > 0, 6„ S K and an infinite variance a-stable random 
variable Yq if and only if the random variable X ~ Xi has a distribution with 
regularly varying tails with index —a G (—2,0), i.e., there exist constants 
p, g with p + q = \ and a slowly varying function L such that 

P(X>x)~p^ and ¥{X -x) ^ , x ^ oo . (1.2) 

This relation is often referred to as tail balance condition. It will be convenient 
to refer to X and its distribution as regularly varying with index a. 

The limit relation (|l.ip is a benchmark result for weakly dependent sta- 
tionary sequences with regularly varying marginal distribution. However, in 
the presence of dependence, conditions for the convergence of the partial sums 
towards a stable limit are in general difficult to obtain, unless some special struc- 
ture is assumed. Early on, a-stable limit theory has been established for the 
partial sums of linear processes {Xt) with iid regularly varying noise with index 
a £ (0,2). Then the linear process [Xt) has regularly varying marginals, each 
partial sum 5'd, d ^ 1, is regularly varying with index a and (Sn) satisfies ()l.ip 
for suitable (a„) and (5„). These results, the corresponding limit theory for the 
partial sums Sn and the sample autocovariance function of linear processes were 
proved in a series of papers by Davis and Resnick (18; 1^ 13) ■ They exploited 



the relations between regular variation and the weak convergence of the point 
processes Nn = X]"=i^a^^Xt' "^hcre £x denotes Dirac measure at x. Start- 
ing from the convergence Nn N, they used a continuous mapping argument 
acting on the points of the processes Nn and N in conjunction with the series 
representation of infinite variance stable random variables. Their proofs heav- 
ily depend on the linear dependence structure. A different, not point process 
oriented, approach was chosen by Phillips and Solo (ji^) who decomposed the 
partial sums of the linear process into an iid sum part and a negligible remain- 
der term. Then the limit theory for the partial sums follows from the one for 
iid sequences with regularly varying marginal distribution. The first result on 
stable limits for stationary processes more general than linear models, assuming 
suitable conditions for non-Gaussian limits, was proved by Davis (jlli). Davis's 
ideas were further developed for mixing sequences by Denkcr and Jakubowski 
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(|22( ) and Jakubowski and Kobus (|34l). The latter paper provides a formula for 
the stable limit for sums of stationary sequences which are m-dependent and 
admit local clusters of big values. A paper by Dabrowski and Jakubowski (fiol ) 
opened yet another direction of studies: stable limits for associated sequences. 

Results for special non-linear time series models, exploiting the structure 
of the model, were proved later on. Davis and Resnick (|2l( ) and Basrak et al. 
Q) studied the sample autocovariances of bilinear processes with heavy-tailed 
and light-tailed noise, respectively. Mikosch and Straumann (|43 ) proved limit 
results for sums of stationary martingale differences of the form Xt = GtZt, 
where {Zt) is an iid sequence with regularly varying ZtS with index a G (0, 2), 
(Gt) is adapted to the filtration generated by (^s)s<t and E|Gt|"+* < oo for 
some (5 > 0. Stable limit theory for the sample autocovariances of solutions 
to stochastic recurrence equations, GARCH processes and stochastic volatility 
models was considered in Davis and Mikosch (13; 14), Mikosch and Starica (|43h . 
Basrak et al. (0); see the survey papers Davis and Mikosch (15; litllTh. 

The last mentioned results are again based on the weak convergence of the 
point processes Nn = StLi^a^^Xt combination with continuous mapping 
arguments. The results make heavy use of the fact that any a-stable random 
variable, a G (0, 2), has a series representation, involving the points of a Poisson 
process. A general asymptotic theory for partial sums of strictly stationary 
processes, exploiting the ideas of point process convergence mentioned above, 
was given in Davis and Hsing (jl2f ). The conditions in Davis and Hsing (12) 
are relatively straightforward to verify for various concrete models. However, 
the a-stable limits are expressed as infinite series of the points of a Poisson 
process. This fact makes it difficult to identify the parameters of the a-stable 
distributions: these parameters are functions of the distribution of the limiting 
point process. 



Jakubowski (13 ll : l33f ) followed an alternative approach based on classical 



blocking and mixing techniques for partial sums of weakly dependent random 
variables. A basic idea of these papers consists of approximating the distribu- 
tion of the sum a~^Sn by the sum of the iid block sums {a~^ Sjni)i=i fc„ such 

that kn = [n/m] — >■ oo and Smi = Sm- Then one can use the full power of 
classical summation theory for row sums of iid triangular arrays. It is also pos- 
sible to keep under control clustering of big values and calculate the parameters 
of the a-stable limit in terms of quantities depending on the finite-dimensional 
distributions of the underlying stationary process. Thus the direct method is in 
some respects advantageous over the point process approach. 



At a first glance, the conditions and results in Jakubowski (j3ll : |33( ) and 



Davis and Hsing 

112) 

look rather different. Therefore wc shortly discuss these 
conditions in Section [2] and argue that they are actually rather close. Our main 
result (Theorem 13. ip is given in Section [3l Using an argument going back to 
Jakubowski (31; .331 ) , wc provide an a-stablc limit theorem for the partial sums 
of weakly dependent infinite variance stationary sequences. The proof only 
depends on the characteristic functions of the converging partial sums. The 
result and its proof are new and give insight into the dependence structure of 
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a heavy-tailed stationary sequence. In Section [321 we discuss the conditions of 
Theorem lS.ll in detail. In particular, we show that our result is easily applicable 
for strongly mixing sequences. In Section 2] we explicitly calculate the parame- 
ters of the Qf-stable limits of the partial sums of the GARCII(1, 1) process and 
its squares, solutions to stochastic recurrence equations, the stochastic volatility 
model and symmetric a-stable processes. 

2 A discussion of the conditions in a-stable limit 
theorems 

2.1 Regular variation conditions 

We explained in Section [T] that regular variation of X with index a e (0,2) in 
the sense of ()1.2p is necessary and sufficient for the limit relation (jl.ip with an 
a-stable limit Ya for an iid sequence {Xt). The necessity of regular variation of 
X with index a G (0,2) in the case of dependent XiS is difficult to establish 
and, in general, incorrect; see Remark 13.21 It is, however, natural to assume 
such a condition as long as one takes the conditions for an iid sequence as a 
benchmark result. 

Davis and Hsing (fl^ assume the stronger condition that the strictly sta- 
tionary sequence {Xt) is regularly varying with index a € (0,2). This means 
that the finite-dimensional distributions of {Xt) have a jointly regularly varying 
distribution in the following sense. For every d ^ 1, there exists a non-null 

Radon measure fid on the Borcl cr-field of M \{0} (this means that /i^ is finite 
on sets bounded away from zero). M = R U {±C!o}, such that 

nP(a-i(Xi,...,Xd)e.)^Md(-), (2.1) 

where A denotes vague convergence (see Kallenberg 
(a„) satisfies 

nP(|X| > a„) - 1. (2.2) 

The limiting measure has the property pLd{xA) = x^'^ ^d{A), t > Q, for Borel 
sets A. We refer to a as the index of regular variation of {Xt) and its finite- 
dimensional distributions. Note that Theorem 3 in (|32h provides conditions 
under which regular variation of the one-dimensional marginals implies joint 
regular variatio n (12.11). 

Jakubowski (j3ll : 1331 ) does not directly assume regular variation of X . How- 
ever, his condition Ul requires that the normalizing sequence (a„) in (|l.ip is 
regularly varying with index 1/a. In ((ssh he also requires the conditions 'T+{d) 
and T_(o?), c? ^ 1, i.e., the existence of the limits 

\im nV{Sd> On) ^h+{d) and Van nV{Sd ^ -a^) ^ b^{d) , d^l. 

n^oo 71^00 

If b+{d) + b^{d) > 0, the regular variation of (a„) with index 1/a is equivalent 
to regular variation of Sd with index a; see Bingham et al. (|3). Condition U2 in 



33), Resnick Um) and 
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(|33[ ) restricts the class of all regularly varying distributions to a subclass. The 
proof of Theorem 13.11 below shows that this condition can be avoided. 

Remark 2.1. Condition (|2.3p is automatically satisfied for regularly varying 
{Xt), where 

b±id) = Md({x e S'' : ±(a;i + ■ ■ ■ + Xd) > I}) . (2.4) 

Since /.i^ is non-null for every d ^ 1 and fj.d{tA) ~ t^°'fi{A), t > 0, we have 
b-f.{d) + b^{d) > 0, d ^ 1. Since (a„) is regularly varying with index 1/a it 
then follows that Sd is regularly varying with index a for every d ^ 1. Since 
(a„) satisfies relation (|2.2[) it then follows that 5+(l) = p and 6-(l) = q with p 
and g defined in equation (jl.2p . In particular p + (7 = 1. The coefficients 6+((i) 
and 6_((i) for d > 1 can be considered as a measure of extremal dependence 
in the sequence {Xt). The two benchmarks are the iid case, b+{d) — pd and 
b-{d) = qd and the case Xi ~ X for all i, &+(d) = pd"' and &_(d) = qd". 

Regular variation of a stationary sequence (Xt ) is a well accepted concept in 
applied probability theory. One of the reasons for this fact is that some of the 
important time scries models (ARMA with regularly varying noise, GARCH, 
solutions to stochastic recurrence equations, stochastic volatility models with 
regularly varying noise) have this property. Basrak and Sogers (@) give some 
enlightening results about the structure of regularly varying sequences. In what 
follows, we will always assume: 

Condition (RV): The strictly stationary sequence (Xt) is regularly varying with 
index a E (0, 2) in the sense of condition (|2.ip with non-null Radon measures 
Hd, d^ 1, and (a„) chosen in (|2.2p . 



2.2 Mixing conditions 



Assuming condition (RV), Davis and Hsing (jl2l ) require the mixing condi- 
tion A{an) defined in the following way. Consider the point process Nn = 
Y^t^i £Xt/a„ s-nd assume that there exists a sequence m = m„ — >■ 00 such that 
kn = [n/rrin] — > 00, where [x] denotes the integer part of x. The condition 
A{an) requires that 

Ec-/-^''^"- (Ec-//''^-)'" ^0, (2.5) 

where / belongs to a sufficiently rich class of non-negative measurable functions 
on R such that the convergence of the Laplace functional Ec " ^ for all / 
from this class ensures weak convergence of (Nn). Relation (|2.5p ensures that 
N„ can be approximated in law by a sum of fc„ iid copies of Nm , hence the weak 
limits of (Nn) must be infinitely divisible point processes. 

The condition ^(a„) is difficult to be checked directly, but it follows from 
standard mixing conditions such as strong mixing with a suitable rate. For 
future use, recall that the stationary sequence (Xt) is strongly mixing with rate 
function (ah) if 

sup \F{An B) -P{A)F{B)\ = 0, h ^ 00 . 

Aedr(...,x_i,Xo),Becr(Xh,Xh+i,...) 
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Jakubowski (j3l[ ) showed that (jl.ip with &„ = and regularly varying (a„) 



implies the condition 



max 

l<k,l<n ,k+l<r 



— >■ , n — >■ oo, X G 



which is satisfied for strongly mixing (Xt). We also refer to the discussion in 



Sections 4-6 of (|33f ) for alternative ways of verifying (|2.6p . Under assumptions 
on the distribution of X more restrictive than regular variation it is shown that 
(jl.ip implies the existence of a sequence Z„ — > oo such that for any kn = o(ln) 
the following relation holds 

/'Ee"'="''°«'^"))''" -Ee"""'^" ^ 0, xeR. (2.7) 



It is similar to condition (|2.5p at the level of partial sums. 

We will assume a similar mixing condition in terms of the characteristic 
functions of the partial sums of (Xt). Write 

ipnj{x) =Ee"''"''^^ , j^l,2,..., (^„ = (p„„, xeR. 

Condition (MX). Assume that there exist m = rUn — > oo such that fc„ = 
[n/m] — > and 

\(Pn{x) - {<fnm{x))''"-\^Q, 71 ^ OO , X E R. (2.8) 

This condition is satisfied for a strongly mixing sequence provided the rate 
function {ah) decays sufficiently fast; see Section [3.2.41 But p.Sp is satisfied for 
classes of stationary processes much wider than strongly mixing ones. Condi- 
tion (MX) is analogous to ^(a„). The latter condition is formulated in terms 
of the Laplace functionals of the underlying point processes. It is motivated by 
applications in extreme value theory, where the weak convergence of the point 
processes is crucial for proving limit results of the maxima and order statis- 
tics of the samples Xi,. . . ,X„. Condition (MX) implies that the partial sum 
processes {a~^Sn) and (a~^ X]f=i ^mi) have the same weak limits, where Smi, 
i = 1, . . . , are iid copies of S,n- This observation opens the door to clas- 
sical limit theory for partial sums based on triangular arrays of independent 
random variables. Since we are dealing with the limit theory for the partial 
sum process {a~^Sn) condition (MX) is more natural than .4(a„) which is only 
indirectly (via a non-trivial continuous mapping argument acting on converging 
point processes) responsible for the convergence of the normalized partial sum 
process (a~^S'„). 

2.3 Anti-clustering conditions 

Assuming condition (RV), Davis and Hsing (flil ) require the anti-clustering 
condition 



lim hmsupPI max \Xi\>xan\\Xo\>xanj=0, x>0, (2.9) 
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where, as before, m = m„ — >■ c« is the block size used in the definition of the 
mixing condition ^(a„). It follows from recent work by Basrak and Segers (@) 
that the index set {« : d < |«| < m„} can be replaced by {z : d ^ i ^ m-n}, 
reducing the efforts for verifying (j2.9|) . With this modification, a sufficient 
condition for (|2.9p is then given by 

lim limsupn y^P(|Xi| > a;a„ , |Xo| > .Ta„) = , x>Q. (2.10) 

I— a 

Relation (|2.10|) is close to the anti-clustering condition D'{xan) used in extreme 
value theory; sec Lcadbcttcr ct al. Leadbetter and Rootzen Jsoh and 

Embrechts et al. (1241 ) . Chapter 5. 

An alternative anti-clustering condition is (38) in Jakubowski (j33l ): 

n-1 

lim lim sup lim sup ^^("- ~ /i)IP(|-^o| > x Un , \Xh\ > x an) = . (2.11) 

h—d 

Assuming regular variation of X and defining (a„) as in (|2.2p . we see that (|2.1ip 
is implied by the condition 



n-1 

lim lim sup lim sup n N P(|X;i| > a;a„ , IXqI > a;a„) = , 



x— >oo n— foo , , 
h—d 



which is close to condition (j2.10p . 

For our results we will need an anti-clustering condition as well. It is hidden 
in assumption (jACp in Theorem 13. II sec the discussion in Section r3. 2.31 

2.4 Vanishing small values conditions 

Davis, Hsing, and Jakubowski prove convergence of the normalized partial sums 
by showing that the limiting distribution is infinitely divisible with a Levy triplet 
corresponding to an a-stable distribution. In particular, they need conditions 
to ensure that the sum of the small values (summands) in the sum a~^Sn does 
not contribute to the limit. Such a condition for a dependent sequence (Xt) is 
often easily established for a e (0, 1), whereas the case a £ [1, 2) requires some 
extra work. 

Davis and Hsing (jl2i ) assume the condition (3.2): 



lim lim sup I 

e ^0 JO ^(-yri 



> X a„J =0, a; > , 
(2.12) 



for a € (0, 2). For an iid sequence (Xt) the relation 

median(a^^S'„) - a^^nE[X/{|x|^ea„}] -> , e > , 

holds. Therefore a~^nKX I^\x\^ea„} ^^re the natural centering constants for 
a~^Sn in stable limit theory. In the case of dependent (Xt), the choice of the 
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latter centering constants is less straightforward; it is dictated by truncation of 
the points in the underlying weakly converging point processes. 



The analogous condition (35) in Jakubowski (j33f ) reads as follows: for each 
a; > 

n 

limlimsuplimsup/''P( ^[Xi "EX/{|x|^e/a^}] >xlan] =0. 

~ ^ ^ (2.13) 

As shown in (1331 ) . this condition is automatically satisfied for a £ (0, 1). 

In our approach the anti-clustering condition (AC) (see below) is imposed 
on the sum of "small" and "moderate" values and we do not need to verify 
conditions such as (|2.12p and (|2.13p . 

3 Main result 

In this section we formulate and prove our main result. Recall the regular 
variation condition (RV) and the mixing condition (MX) from Sections 12.11 
and 12.21 We will use the following notation for any random variable Y: 

Y= (yA2) V(-2). 

Notice that |F| = |r| A 2 is subadditive. 

Theorem 3.1. Assume that (Xt) is a strictly stationary process satisfying the 
following conditions. 

1. The regular variation condition (RV) holds for some a G (0,2). 

2. The mixing condition (MX) holds. 

3. The anti- clustering condition 



limlimsup— V ¥.\x On^iS, - So) x aV,^ Xi =0, 2;eM, (AC) 
d->oo „_>.oo m .-^^ I 

holds, where m = ?7!,„ is the same as in (MX). 
4-. The limits 

lim (5+(rf)-&+(d- 1)) = c+ and lim (6_ (d) - 6_ (d - 1)) = c_ , (TB) 

exist. Here hji-[d),b^{d) are the tail balance parameters given in (|2.3p . 
5. For a > 1 assume E(Xi) = and for a = 1, 

lim limsupn|E(sin(a^^S'<j))| = 0. (CT) 

d—^oo fi — 
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Then c+ and c_ are non-negative and {a~^Sn) converges in distribution to an 
a-stable random variable (possibly zero) with characteristic function ipa{x) = 
exp(— c+, c_)), where for a ^ 1 the function Xa{x,c+,C-),x G R, is 
given by the formula 

— — f (c+ + c_) cos(7ra/2) — i sign(a;)(c+ — c_) sin(7r a/2)) , 
1 — a V / 

while for a ~ 1 one has 

Xi(x, c+, c_) ~ 0.5 7r(c+ + c_) + i sigii(a;) (c+ — c_) log x E M. 
We discuss the conditions of Theorem 13.11 in Section 13.21 In particular, 



we compare them with the conditions in Jakubowski pit l33f ) and Davis and 



Hsing (|l2r). If the sequence (X„) is mo-dependent for some integer mg ^ 
1, i.e., the cr-fields a{. . . ,X-i,Xq) and a{X„i„+i, X„ig^2, ■ ■ •) are independent, 
the conditions (MX), (AC) and (TB) of Theorem 13.11 are automatic; see 
Section 14.11 The surprising fact that c+ and c_ arc non-negative is explained 
at the end of Section [3.2.21 

Remark 3.2. Although Theorem 13.11 covers a wide range of strictly stationary 
sequences (see in particular Section 2]) condition (RV) limits the applications 
to infinite variance (in particular unbounded) random variables X„. The referee 
of this paper pointed out the surprising fact that there exist strictly stationary 
Markov chains (Yn), suitable bounded functions / and a sequence (a„) with 
a„ = n^^°'£{n) for some slowly varying function £ such that the sequence of 
the normalized partial sums (a^^S'„) of the sequence (A„) = {f{Yn)) converges 
in distribution to an infinite variance stable random variable. Then (RV) is 
obviously violated. Such an example is contained in Gouezel ([29|). Theorem 1.3. 

Other examples of stable limits for sums of bounded stationary random vari- 
ables (of different nature - non-Markov and involving long-range dependence) 
are given in (fs^). Theorems 2.1 (ii) and 2.2 (ii). 

Remark 3.3. It might be instructive to realize that in limit theorems for weakly 
dependent sequences properties of finite dimensional distributions can be as bad 
as possible. For example it is very easy to build a 1-dependent sequence having 
no moment of any order and such that its (centered and normalized) partial 
sums still converge to a stable law of order a € (0, 2]. Let 

Xn — -\- Sn ^n—li 

where {Yn) is an iid sequence of a-stable random variables and (e„) is an iid 
sequence without any moments (that is £^(|£o|") = +oo for any positive a) and 
the two sequences are independent. Then the (centered and normalized) partial 
sums of (A„)) have the same limit behavior (in distribution) as that of (Yn). 

3.1 Proof of Theorem ISTTl 

For any strictly stationary sequence (Xt)t£Z it will be convenient to write 
^0 = 0, Sn^Xi + --- + Xn, = A_„ + --- + A_i, n^l. 
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Let Smi, i = 1, 2, . . . , be iid copies of Sm- In view of (|2.8p the theorem is proved 
if we can show that {a~^ 'Y^^Zi Smi) has an a-stable limit with characteristic 
function ipa. For such a triangular array, it is implied by the relation 

kn{ipn7n{x) - I) ^ \0gija{x) , X&R. (3.1) 

Indeed, notice first that the triangular array {a~^ Smi)i=i,...,k„ of iid random 
variables satisfies the infinite smallness condition. Then apply Lemma 3.5 in 
Petrov (|46[ ) saying that for all 2; G M and sufficiently large n, 

l0g<^,i,„(a;) ^ (Pn,n{x) - 1 + On,n{(p„,n{x) - 1)^ , 

where |6'„,„| ^ 1. Thus 

|fcn(Vrim(a;) - 1) - fc„ log i^„,„ (a;) | ^ fc„ \(pn,nix) - 1^ ^ cfc,7^ ^ . 

Here and in what follows, c denotes any positive constants. Our next goal is to 
find a suitable approximation to the left-hand side in p.ip . 

Lemma 3.4. Under (RV) and (AC) the following relation holds: 

kniVnmix) - I) - n {(pnd{x) - (Pn,d-l{x)) =0, XER. (3.2) 

Moreover, if (Xn) is mo -dependent for some integer mo ^ 1, then 



lim limsup 



lim 

n— J-cxD 



kn{(pnmix)~l)-n{ipnd{x)-ipn,d-i{x)) =0, xGR, d > mQ . (3.3) 



The proof is given at the end of this section. By virtue of (RV) and (|2.3p . Sd 
is regularly varying with index a £ (0, 2); see Remark |2. II Therefore it belongs 
to the domain of attraction of an a-stable law. Theorem 3 in Section XVII. 5 of 
Feller (|26h yields that for every d ^ 1 there exists an a-stable random variable 
Za{d) such that 

an^y!('^d* - e„d) ^ ^a(c^) where e„d = < ' (3.4) 

^ [E(sm(5d/a„)) a = I. 

The limiting variable Za{d) has the characteristic function 

V'a,d(a;) = c-x:p{-\x\°'xa{x,h+{d),h^{d))) , x e R . 



Applying Theorem 1 in Section XVII. 5 of Feller ()260 . we find the equivalent 
relation 

n(v3„d(a;)e-"^"^^ - 1) Xogi^^A^) ^ xgM, 
and exploiting condition (jTBp . for a; e R, 

n(<p„d(a;)e-^^"''^ - ¥'«,d-i(a;)e ""'"■''-i^") 

-> \og%l^aA{x) - \og%j}aA-i{x) as n cxo (3.5) 

— > log V'q (x) as (i — > oo . 
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For a 7^ 1 we have end = 0. Therefore p.2p impUes 

This finishes the proof in this case. For a = 1 we use the same arguments but 
we have to take into account that end does not necessarily vanish. However, we 
have 

\Vnd{x)-^nd{x)e-'^"'^\ =^ |l-e-*^"<^n 
and using (jCT[) . we obtain 

lim limsupn |(<^„d(a;) - fn,d-^i{x)) - {'^nd{^)e~"'"'''^ - <(9„,d_i (2;)e ""'"•''-i"'') | = 0. 

This proves the theorem. dProof of Lemma \3.4\ Consider the following 

telescoping sum for any n, m n and d < m: 

m—d 

Ifinmix) - 1 = (Pnd{x) " 1 + (V'n,d+j{x) - (pn,d-l+j {x)) . 

J = l 

By stationarity of (Xt) we also have 

m {ipnd{x) - ipn,d-i{x)) = 

m — d 
j = l 

Taking the difference between the previous two identities, we obtain 

(V9„m(x-) - 1) - m {<find{x) - (p„^d-lix)) = 

-(d- 1) - 1) + d{ifn.d^i{x) - 1) + 

m — d 

By stationarity, for any fc ^ 1, fnk{x) = Ee*^""^^-'=. Therefore any summand 
in the latter sum can be written in the following form 

= E^e"""^'^"''"^ (l — e"'^""^"^"^) (l — e ~*'^'^"^'^"''"^)^ . 

Using the fact that x — > exp(ix) is a 1-Lipschitz function bounded by 1, the 
absolute value of the expression on the right-hand side is bounded by 

m\xa-^S.,\A2){\xa-'X^d~,\h2)). 



11 



Collecting the above identities and bounds, we finally arrive at the inequality 

\k„ {ipnm(x) - 1) - n {(pnd{x) - <^„,d_ i (z)) | 

^ kn (d- 1) \(p„d{x) - 1| + fc„ d |<^„,d-i(a;) - 1| 



+kn ^ E |a;a„"^(S'j - Sd) xun^Xi . 

j=d+l 

The last term on the right-hand side converges to zero in view of assumption 
(|ACp when first ?i — > oo and then d — >■ oo. In the mg-dependent case, the last 
term on the right-hand side converges to zero whenever d > toq and n — >■ oo. 
To prove that the first two terms also converge to zero, let us notice that, under 
(RV), n{Lpnd{x) - 1) ^ Xa{x,b+{d),b^{d)) and so 

lim kn{d - l){ipnd{x) - 1) = lim TO„~^(d - l)n((^„rf(x) - 1) = 0. 

n— J-oo n— foo 

This proves the lemma. □ 

Remark 3.5. Balan and Louhichi (Q) have taken a similar approach to prove 
limit theorems for triangular arrays of stationary sequences with infinitely di- 
visible limits. Their paper combines ideas from Jakubowski (|33l ). in particular 
condition (jTBp . and the point process approach in Davis and Hsing ()l2l ). They 
work under a mixing condition close to A{an). One of their key results (The- 
orem 2.6) is the analog of Lemma [3.41 above. It is formulated in terms of the 
Laplace functionals of point processes instead of the characteristic functions of 
the partial sums. Then they sum the points in the converging point processes 
and in the limiting point process to get an infinitely divisible limit. The sum 
of the points of the limiting process represent an infinitely divisible random 
variable by virtue of the Lcvy-Ito representation. As in Davis and Hsing (fT^ 
the method of proof is indirect, i.e., one does not directly deal with the partial 
sums, and therefore the results are less explicit. 

3.2 A discussion of the conditions of Theorem 13.11 

3.2.1 Condition 5 

It is a natural centering condition for the normalized partial sums in the cases 
a = 1 and a G (1, 2). In the latter case, E\X\ < oo, and therefore EX = can 
be assumed without loss of generality. As usual in stable limit theory, the case 
a = 1 is special and therefore we need condition (jCTp . It is satisfied if Sd is 
symmetric for every d. 

3.2.2 Condition ((TB|) 

If C-|_ + c_ =0 the limiting stable random variable is zero. For example, assume 
Xn =: Yn — Yn-i for an iid regularly varying sequence (y„) with index a € (0, 2). 
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Then Sd = Yd — Yq is symmetric and regularly varying with index a. By 
the definition of (a„), = b-{d) = 0.5, hence c_|- = c_ = 0. Of course, 

a-'^Sn 4 0. 

In the context of Theorem 3.1 in Jakubowski (fssl ) (although the conditions 
of that result are more restrictive as regards the tail of X) it is shown that 
(|TBp is necessary for convergence of (a~^Sn) towards a stable limit. Condition 
(|TB[) can be verified for various standard time series models; see Section 21 
The meaning of this condition is manifested in Lemma 13.41 It provides the link 
between the regular variation of the random variables Sd for every d ^ 1 (this 
is a property of the finite-dimensional distributions of the partial sum process 
{Sd)) and the Levy measure I'a of the a-stable limit. Indeed, notice that (|TB[) 
implies that, for every x > 0, with fe+(0) = 0, 

= lim ;.5](&+(i)-5+(*-l))x-« 

2—1 

1 

= lim —/y lim n (P(S'i > a; a„) — P(S'i_i > a; a„)) , 

and a similar relation applies to c_a::~". Then 

Va(x^(X}) = Ca^x"'^ and i^Q(— oo, — x) = c_ a::^" , x>0, 

determine the Levy measure of the a-stable limit distribution with the char- 
acteristic function i/iq, given in Theorem 13.11 In particular. Lemma 13.41 implies 
that as n — > oo 

P(S'm > X Qn) ^ Vai^X, Oo), a; > f3 6) 

fc„ P(5'm -Xa„) ^- Va{-(X),-x\, 

The latter relation opens the door to the limit theory for partial sums of trian- 
gular arrays of nd copies (S'mi)i=i,...,fc„, of Sm- Notice that the relations p.6p 



are of large deviations type in the sense of (|33r). We refer to (|35f ) for their 
multi-dimensional counterparts. 

Let us notice that although one cannot ensure that h±{d) ^ h±{d — 1) ^ 
for sufficiently large d, the constants c+,c- are non- negative. It is immediate 
from the observation that 



c± = lim (6±(d) -6±(d- 1)) 

d—foo a ^ — ' d^-oo a 
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Remark 3.6. Recall the two benchmark examples of Remark l2.1l If (Xt) is an 
iid sequence regularly varying with index a > the limits ~ p and ~ q 
always exist and conditions (MX), (AC) are automatically satisfied. Then, 
under (CT), we recover the classical limit results for partial sums with a-stable 
limit. On the other hand, ii Xi = X for all i, then c+ = c_ = if < a < 1, 
c+ ~ p and c_ = q ii a = 1 and c+ and c_ are not defined otherwise. This 
observation is in agreement with the fact that a~^Sn = n^~^^"£{n)X for some 
slowly varying function £. 



3.2.3 Sufficient conditions for ((AC|) 



Condition (jAC[) is close to the anti-clustering conditions in (jl2t 1331 ) discussed 
in Section [2.31 In what follows, we give some sufficient conditions for (|AC[) . 
These conditions arc often simple to verify. 

Lemma 3.7. Assume the conditions of Theorem \3. 1\ and that (Xt) is strongly 
mixing with rate function (ah)- Moreover, assume that there exists a sequence 
r„ — >■ oo such that r„/m„ — > 0, nar„ — > and one of the following three 
conditions is satisfied. 



or 



or 



limd_+oohmsup„^^n [J2i=d+i^i\^i\ >an,\Xi \ > a„) 

+ P{\j:Zd+iX,J{\x^l^aJ>an,\Xi\>a„)] - 0. 



lim limsupnP( max \Xi\ > a„/r„ , \Xi\ > a„) = . 

d— >oo n— >oo i—d+l,...,rn 



lim limsupnP(|S'r„ — Sd\ > fln , |-^i| > Qn) = . 



(3.7) 

(3.8) 
(3.9) 



Then ((AC|) holds. 



Proof. Let us recall that the function y i-^ \y\ is subadditive. We decompose 
the sum in ()AC[) as follows. 



\ J2 + \^ xOn'^iSj - Sd) xan^Xi = Ji{n) + J2{n) . 

\]=d+l J=r„ + l/ 

We will deal with the two terms Ji{n) and J2{n) in different ways. For the sake 
of simplicity we assume a; = 1 . 
We start by bounding J2{n). 



J2{n) ^ fc„ ^ [e 

i=i-„+i 



0-n^{Sj ~ Sr,^) On^ Xi + an^{Sr^—Sd) 



j=r„+l 



nE 



i^{^r„ — Sd) On^Xi 



■hlin) + J22{n) . 
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We bound a typical summand in J2i{n), using the strong mixing property 

E 



^ car 



Moreover, we have for j > r„ 



hence 



E 



sCj E 



O-n^ Sj-r„ + l E Un^Xi 

Thus we arrive at the bound 

J2ii'n') ^ cnar„ +cnm (m 



ttn^Xi 



a„ ^Xi 



Observe that 



E\an'X,\ = E(a,-i|Xi|/^^-X|^^|^2jA2/^^-.|^^|>2}) 

< 2P(|Xi| >2a„)- 
Therefore, by definition of (a„) and since na^,, — >■ by assumption, 
J2i(7i) = 0{narj + 0{m/n) = o(l) . 

We also have 

J22(n) =^ cnF{\Sr„- Sd\>2an,\Xi\>2an) 

^ cn 2 P(|X,| > a„,|Xi| > 2a„) 



> 2a„ , |Xi| > 2a„ 



and 



J22N cP( max > 2 a„/r„ I > 2 a„). 

i—d-yX^.-.^rn 



Thus, under any of the assumptions (|3.7p - (|3.9p . hmf;_^oo l™sup^j_^^ J^iij^ — 0. 
Finally, 

Ji{n) ^ cK £ P(|^, - Sd\ > 2a„ , > 2a„) 
j=d+i 

<^ !^nP(|Xi| > 2a„) = o(l). 
m 
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Collecting the bounds above we proved that (|ACp holds. □ 
3.2.4 Condition (MX) 

As wc have already discussed in Sections 12.11 and 12. 2[ the condition (MX) is 
a natural one in the context of stable limit theory for dependent stationary 
sequences. Modifications of these conditions appear in Davis and Hsing (|l2l ). 
Jakubowski (31; 33|). We also discuss the existence of sequencees m = m„ — )■ oo 
and r = r„ — > oo such that r/m — > and m/n — !■ to be used in Lemma 13.71 in 
the context of strong mixing. 

Lemma 3.8. Assume that (Xt) is strongly mixing with rate function (ah). In 
addition, assume that there exists a sequence e„ — >■ satisfying 

""[^^(ayriAr,)] ^ . (3.10) 

Then (MX) holds for some m ~ m„ —> oo with fc„ ~ [n/m] — > oo. Moreover, 
writing rn = [en(art A n)], then r„/m„ — s- for this choice of (to„) and 

nar^-^O. (3.11) 

If the tail index a ^ 1, then p.ll|) turns into na^^^^n] ~^ which is more 
restrictive than p.lOp . On the other hand, if the tail index a is close to 2, (|3.10p 
is not implied by polynomial decay of the coefficients ah- Then a subexponential 
decay condition of the type a„ ^ C exp(— cn^) for some C, c, 6 > implies p.lOp . 
and then ()3.1ip follows. 

Proof. We start by showing that (|2.8p holds for a suitable sequence (m„). Let 
fnmS be the characteristic function of a^^ X^i^i Um~s.i for some 5 = (5„ and 

Uji = ^ Xk 

k={i-l)j+l 

a block sum of size j. Using that characteristic functions are Lipschitz functions 
bounded by 1 and writing q ~ fc„m, for .t e M, 



\y:>q{x) - (pnm5{x)\ ^ E(\ ^ [^5, 



k„ 

X 



j = l 



A2 



^ E(^|Ai|A2 



/ Y{\Xi\> man/{5n)s)ds. 



The right-hand side approaches zero if mnan/{Snn) ^ oo as n — > oo. Under 
this condition, the same arguments yield 

li'Pnmix))''" - {(p„,m~5{x))''"-\ ^ and | (^g (.t) - (/3„ (x) | -> , 
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as soon as anjm — )■ oo. Next we use a standard mixing argument to bound 

n — m.mo 

+ \ip„i-s{x)ipn-m,msix) - „_5 (x)) I™/"' | . 

The first term on the right-hand side is the covariance of bounded Lipschitz 
functions of Sm-s and Sn — S„i- Hence it is bounded by as- Iterative use of 
this argument, recursively on distinct blocks, shows that the right-hand side is 
of the order {n/m)as- Thus we proved that (|2.8p is satisfied if 

n/mas^O, man/{Sn)-^oo and a„/?Ti — > oo . (3-12) 

Choose m„ = [^/e^{an A n)], S = [m'^/n] and assume p.lOp . Then p. 121) holds, 
(|2.8p is satisfied and m/n ^ \/S/n — > 0. 

Finally, if p.lOp is satisfied choose r„ = [e„(a„ A n)]. Then nar^ — > and 
rn/m — > are automatic. □ 

4 Examples 

4.1 mo-dependent sequences 

Consider a strictly stationary sequence (Xn) satisfying condition (RV) and 
which is mo-dcpcndcnt for some integer mo ^1. In this case, ah = for 
h > mo- Then, by virtue of Lemma 13.81 condition (MX) is satisfied for any 
choice of sequences (m„) such that m„ — > oo and 7ti„ — o(n). Moreover, (AC) 
follows from Lemma 13.71 for any (r„) such that r,i — > oo and r„ = o(77i„). We 
verify the validity of condition (j3.8p . Then for (r,i) growing sufficiently slowly, 

nP( max \Xi\ > a„/r„ ,\Xi\ > a„) 

s$ nrnP{\Xi\> an)P{\Xi\> an/rn) 
= 0(r„P(|Xi| >a„/r„)) = o(l). 

Thus, in the mo-dependent case we have the following special case of Theo- 
rem [STTJ 

Proposition 4.1. Assume that (Xt) is a strictly stationary mo-dependent se- 
quence for some mo ^ 1 which also satisfies condition (RV) for some a € (0, 2). 
Moreover, assume E(Xi) = for a > 1 and Xi is symmetric for a = 1. 
Then the conclusions of Theorem \3.1\ hold with Cj^ — b-i.(mo -|- 1) — 5_|_(mo) and 
c_ = 5_(m.o + 1) - &-(mo). 

Proof. We have already verified conditions (MX) and (AC) of Theorem 13.11 
Following the lines of the proof of Theorem 13.11 with Cnd — 0, we arrive at p.Sp 
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for every d ^ 1. In view of the second part of Lemma the right-hand side of 
(|3.5p is independent of d for d > rriQ as the hmit of kn{ipnmix) — 1) as n — >■ oo. 
This finishes the proof by taking d = niQ + 1. □ 

Wc mention in passing that we may conclude from the proof of Proposi- 
tion |43] that condition (TB) is satisfied since c+ = b+{d + 1) — b+{d) and 
c_ = b-{d + 1) — b-{d) for d > mo. This is a fact which is not easily seen by 
direct calculation on the tails oi Sd, d > mo- 

4.2 The stochastic volatility model 

The stochastic volatility model is one of the standard econometric models for 
financial returns of the form 

Xt = at Zt , 

where the volatility sequence (at) is strictly stationary independent of the iid 
noise sequence (Zt). See e.g. Andersen et al. (f!]) for a recent reference on 



stochastic volatility models or the collection of papers (|5lh. 



Conditions (RV), (|TB|) and (|CT)) 

We assume that Z is regularly varying with index a > 0, implying that {Zt) is 
regularly varying. We also assume that Ea^ < oo for some p > a. Under these 



assumptions it is known (see Davis and Mikosch (|14I)) that (Xt) is regularly 



varying with index a, and the limit measure in (|2.1I) is given by 

d 

^jid{dxi,...,dxd) = '^\a{dxi)'^eo{dx.j) , (4.1) 

i=l i^j 

where is Dirac measure at x, 

\a{x,(X)) = px~" and Aq(— oo, — a;] = , x>0, 

and 

n = luTL „, and q = Imi ^, , 4.2 

^ x^o^f{\Z\>x) ^ x^oo P(|Z| > X) ' ^ ' 

are the tail balance parameters of Z. This means that the measures /id are 
supported on the axes as if the sequence {Xt) were iid regularly varying with 
tail balance parameters p and q. By virtue of (|4.ip and (|2.4p we have 6+(rf) = pd 
and 6_ (d) = g d, hence c+ = p and c_ = q. We also assume EZ = for a > 1 . 
Then EX = 0. If a = 1 we assume Z symmetric. Then Sd is symmetric for 
every d ^ 1 and (jCXp is satisfied. 
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Conditions (MX) and ((AC|) 



In order to meet (MX) we assume that (oj) is strongly mixing with rate function 
(ah)- It is well known (e.g. Doukhan (|23f)) that (Xf) is then strongly mixing 
with rate function (4a/i). 

It is common use in financial econometrics to assume that (log at) is a Gaus- 
sian linear process. The mixing rates for Gaussian linear processes are well 
studied. For example, if (log at) is a Gaussian ARM A process then (ah) decays 
exponentially fast. We will assume this condition in the sequel. Then we may 
apply Lemma 13.81 with r„ = n'^'^ , m„ = rf'^ , < 71 < 72 < 1 for sufficiently 
small 71 and 72, to conclude that (MX) holds and nar„ 0. 

Next we verify (jACp . We have by Markov's inequality for small e > 0, 

nP( max \Xi\ > a„/r„, \Xi\ > a„) 

2— c?+l,...,r„ 

< n 2 F{\X,\ > a,Jr„, \Xi\ > a„) 

i=d+l 

< n ^ P(max(crj, CTi) min(|Z,|, |Zi|) > a„/r„) 

< n{r,Ja,,r+^ £ E(max«+^ 

The right-hand side converges to zero if we choose 71 and e sufficiently small. 
This proves p.Sp and by Lemma [5771 also (jACp . 

Proposition 4.2. Assume that (Xt) is a stochastic volatility model satisfying 
the following additional conditions: 

(a) {Zt) is iid regularly varying with index a € (0, 2) and tail balance param- 
eters p and q. 

(b) For a £ (1, 2), KZ = 0, and for a = 1, Z is symmetric. 

(c) (logCTt) is a Gaussian ARMA process. 

Then the stochastic volatility process (Xt) satisfies the conditions of Theorem \3.1\ 
with parameters c^ = p and c^ — q defined in (j4.2p . 

Hence a stochastic volatility model with Gaussian ARMA log-volatility se- 
quence satisfies the same stable limit relation as an iid regularly varying sequence 
with index a € (0, 2) and tail balance parameters p and q. 

In applications it is common to study powers of the absolute values, (|Xt|^), 
most often for p = 1,2. We assume the conditions of Proposition 14.21 Then 
the sequence (X^) is again a stochastic volatility process which is regularly 
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varying with index q;/2 £ (0, 1). It is not difficult to see that the conditions of 
Proposition 14.21 arc satisfied for this sequence with 6_ (d) = and b+ (d) = d, 
hence c+ = 1 and c_ = 0. 

A similar remark applies to with one exception: the centering condi- 

tion (jCTp cannot be satisfied. This case requires special treatment. However, 
the cases a 1 arc similar. For a < 1. (|^t|) is a stochastic volatility model 
satisfying all conditions of Proposition 14.21 For a G (1,2) wc observe that 



a, 

1=1 t=i 



n 

^Y,at{\Zt\-E\Z\) + op{l). 



t=l 



In the last step we applied the central limit theorem to {(Jt)- Then the process 
{(Jt{\Zt\ — E\Z\)) is a stochastic volatility model satisfying the conditions of 
Proposition 14.21 with c+ = 1 and c_ = 0. 



4.3 Solutions to stochastic recurrence equations 

We consider the stochastic recurrence equation 

Xt^AtXt-i+Bt, teZ, (4.3) 

where {{At, Bt)) constitutes an iid sequence of no n- negative random variables At 
and Bt- Various econometric time scries models {Xt) have this form, including 
the squared ARCH(l) process and the volatility sequence of a GARCH(1,1) 
process; sec Section [4.41 The conditions Elog^ < and E|logi?| < oo are 
sufficient for the existence of a strictly stationary causal solution {XA to (|4.3p 
such that {Xn)n^o and {{An, Bn))n^i ai'c independent; see Kesten (|38l). 



Condition (RV) 

Kesten ((ssh and Goldie ((28h showed under general conditions that X has almost 
precise power law tail in the sense that 

P{X > x) ^ Cox-" (4.4) 

for some constant cq > 0, where the value a is given by the unique positive 
solution to the equation 

EA'^ = 1 , K > . 
Wc quote Theorem 4.1 in Goldie ((28h to get the exact conditions for (|4.4p . 

Theorem 4.3. Assume that A is a non-negative random variable such that the 
conditional law of A given A ^ is non- arithmetic and there exists a > such 
that = 1, E(A" log+ A) < oo. Then -oo s$ Elog A < and E{A°' log A) e 
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(0,cx)). Moreover, if KB" < oo, then a unique strictly stationary causal solution 
(Xt) to (|4.3p exists such that (|4.4p holds with constant 

'° = am-' log A) ■ ^'-^^ 

The condition of non-arithmeticity of the distribution of A is satisfied if A has a 
Lebesgue density. In what follows, we assume that the conditions of Theorem\^ 
are satisfied. 

Iterating the defining equation (j4.3p and writing 



= Ai • • • At , t ^ 1 , 

we see that 

= Xo(ni,n2,...,nd) + i?d, (4.6) 



where Rd is independent of Xq. Under the assumptions of Theorem 14.31 the 
moments EA" and EB" are finite, hence E(i?2^) < oo and Vi\Rd\ > x) = 
o(P(|Xn[ > x)). By a muhivariate version of a result of Breiman (|8|) (see Basrak 
et ah (jsl)) it follows that the first term on the right-hand side of (|4.6p inherits 
the regular variation from Xq with index a and by a standard argument (see 



Jessen and Mikosch (|36f). Lemma 3.12) it follows that (Xi, . . . , Xd) and the first 
term on the right-hand side of (|4.6p have the same limit measure ^d- Hence the 
sequence {Xt) is regularly varying with index a, i.e., condition (RV) is satisfied 
for a > with EA" = 1. 



Condition (|TBl) 

Next we want to determine the quantities b+{d). Choose (a„) such that nP{X > 
On) ^ 1, i.e.. On = (cqtt,)^/", and write 

d 

Td = ^ n, , d^i. 

We obtain for every d ^ 1, by (|4.6p . 

nP{Sd > a„) - nP(Xord > a„) n¥{Xo > a„)E(Trf«) ^ E(T^") = b+{d) . 

Here we again used Breiman's result H) for F{XoTd > x) - E(T^")P(Xo > x) 
in a modified form. In general, this result requires that E(r^+*) < oo for some 
S > 0. However, if P{X > x) ~ cox~°', Breiman's result is applicable under 
the weaker condition EA" < oo; see Jessen and Mikosch (fsih. Lemma 4.2(3). 
Of course, b-{d) = 0. We mention that the values b±{d) do not change if Sd is 
centered by a constant. 

Our next goal is to determine c+. Since KA" = 1 we have 

b+{d + 1) - b+{d) = E[(l -t- Tdr ~ T2\ . (4.7) 
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The condition EA" = 1 and convexity of the function = KA'^, k > 0, imply 
that E log A < and therefore 

oo 

r<j Too = Hi < oo . 

i=l 

Therefore the question arises as to whether one may let d — >■ oo in (|4.7p and 
replace in the limit by Too ■ This is indeed possible as the following dominated 
convergence argument shows. 

If a e (0, 1] concavity of the function f{x) = x" yields that (l+rrf)"-T^" ^ 1 
and then Lebesgue dominated convergence applies. If a € (1, 2), the mean value 
theorem yields that 

where ^ g (0, 1). Hence (l + Td)" -T^" is dominated by the fimction a[T^~'^ + l]. 
By convexity of gin), k > 0, we have E{A"~^) < 1 and therefore 

oo oo 

E(T^-i) ^ ^E(nf-i) = ^(E(yl"-i))' =E(^"-i)(l -E(^"-i))-i < oo . 

1=1 i=l 

An application of Lcbcsguc dominated convergence yields for any a G (0, 2) 
that 

c+= lim[6+(d+l)-5+(rf)]=E[(l + roo)"-r^]e(0,oo). (4.8) 

Remark 4.4. The quantity Too has the stationary distribution of the solution 
to the stochastic recurrence equation 

Yt = AtYt-i + l, t € Z . 

This solution satisfies the conditions of Theorem 14.31 and therefore 

P(ro >x)^ P(Too > x) - Ci x^" , 

with constant 

E[yi" - (Ai yo)"] + Ai Yor - (^1 Yor] 



Cl 



a E(A" log A) a E(A" log A) 



In particular, 1L(T^) = oo. This is an interesting observation in view of c+ € 
(0,oo). It is also interesting to observe that the limit relation (|4.8p implies that 

= = E[rf-i/"rd]" ^ E[(l + TooT - T^] , 

although d~^/°'Td 0. This relation yields some information about the rate 
at which Td — >-' T^o- 
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Condition (MX) 

The stationary solution {Xf) to tlie stocliastic recurrence equation (|4.3p is 
strongly mixing with geometric rate provided that some additional conditions 
are satisfied. For example, Basrak et al. (Q), Theorem 2.8, assume that the 
Markov chain (Xt) is /z-irreducible, allowing for the machinery for Feller chains 
with drift conditions as for example explained in Feigin and Tweedie (|25r ) or 
Meyn and Tweedie (0). The drift condition can be verified if one assumes 
that At has polynomial structure; sec Mokkadem (|45[ ). The latter conditions 
can be calculated for GARCH and bilinear processes, assuming some positive 
Lebesgue density for the noise in a neighborhood of the origin; see Basrak et al. 
(0), Straumann and Mikosch (|53h . 

In what follows, we will assume that (Xt) is strongly mixing with geometric 
rate. Then, by Lemma 13. 8[ we may assume that we can choose r„ = n'^'^ , 
m„ = n'''^ for sufficiently small values < 71 < 72 < 1. Then (MX) holds and 
nar„ — > 0. 

Condition (|AC|) 

We verify condition p.9p and apply Lemma [5771 It suffices to bound the quan- 
tities 

In{d)=V{\Sr„-Sd\>an\Xo>a„). 
Writing Ils.t = Yli=s s ^ t and list ~ 1 for s > t, we obtain 

i 

x, = XoU, + Y, n/+i,. Bi = n, + a , 1^1. 

Then, using the independence of Xq and Ci, i ^ 1, applying Markov's inequality 
for K < a A 1 and Karamata's theorem (see Bingham et al. (0)), 

d > an/2) 

i=d+l i=d+l 

n \ ^ i=d+l 1=1 

00 

< c ^ {E{A'')Y + crl+''a-''{\-EA'')-^EB'' 

i=d+l 



< 



c{{nA'')Y + rl+-a--). 



Here we also used the fact that E{A'^) < 1 by convexity of the function g{Hi) = 
E{A'^), K > 0, and g{a) = 1. Choosing r„ = n'''^ for 71 sufficiently small, we see 
that 

lim lim sup /„ (d) = . 
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This proves ([X^ . 



Condition 5 

Since X is non-negative, (jCT[) cannot be satisfied; the case a = 1 needs special 
treatment. We focus on the case a e (1,2). It is not difficult to sec that all 
calculations given above remain valid if we replace Xt by Xt — EX . provided 71 
in r„ = n'^'i is chosen sufficiently small. 
We summarize our results. 

Proposition 4.5. Under the conditions of Theorem \4^.S\ the stochastic recur- 
rence equation (|4.3p has a strictly stationary solution {Xt) which is regularly 
varying with index a > given by E[A"] = 1. If a £ (0,1) U (1,2) and {Xt) is 
strongly mixing with geometric rate the conditions of Theorem \8.1\ are satisfied. 
In particular, 



{cqu) (5„ - 6„) A Za , where 6„ = 



forae (0, 1) , 

ESn=nEX for a e {1,2), 



where the constant cq is given in (j4.5p and the a- stable random variable Z^ has 
characteristic function ipa{t) = exjp{—\t\°'Xa{t, c+ , 0)) , where 

00 

c+ =E[(1 + Too)" -T^] e (0,00) and T^, = ^ ^li • • • A, . 

i=l 

Remark 4.6. Analogs of of Proposition have recently been proved in Guiv- 
arc'h and Le Page (jSOj) in the one-dimensional case and in Buraczewski et al. 
13), Theorem 1.6, also in the multivariate case. The results are formulated for a 
non-stationary version of the process {Xn) starting at some fixed value Xq = x. 
(This detail is not essential for the limit theorem.) The proofs are tailored 
for the situation of stochastic recurrence equations and therefore different from 
those in this paper where the proofs do not depend on some particular structure 
of the underlying stationary sequence. 

4.4 ARCH(l) and GARCH(1,1) processes 

In this section we consider the model 

Xt = at Zt , 

where {Zt) is an iid sequence with EZ = and var(Z) = 1 and 

af = ao + {aiZf_i + f3i) a'i_^ . (4.9) 

We assume that ao > and the non- negative parameters ai. Pi are chosen such 
that a strictly stationary solution to the stochastic recurrence equation (|4.9p 
exists, namely, 

- 00 ^ Elog(aiZ2 +/3i) < 0, (4.10) 
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see Goldie (|28^ . cf. Mikosch and Starica (|43( ). Then the process (Xt) is strictly 
stationary as well. It is called a GARCH(1,1) process if ai/3i > and an 
ARCH(l) process ii /3i ~ and ai > 0. Notice that condition (j4.10p implies 
that /3i e [0, !)• 

As a matter of fact, these classes of processes fit nicely into the class of 
stochastic recurrence equations considered in Section 14.31 Indeed, the squared 
volatility process (af) satisfies the stochastic recurrence equation (|4.3I) with 
Xt = cr^. At = a\Zl_^ + /?! and Bt = a^. Moreover, the squared ARCII(l) 
process (X^) satisfies with Yt = X}, At = aiZf and Bt = ao Zf. 

A combination of the results in Davis and Mikosch (jl3l ) for ARCII(l) and 
in Mikosch and Starica (j43r ) for GARCII(1, 1) with Proposition 23] above yields 
the following. 

Proposition 4.7. Let {Xt) be a strictly stationary GARCH(1,1) process. As- 
sume that Z has a positive density on M and that there exists a > such that 

E[(aiZ2 + = 1 and E[{aiZ'^ + Pi)" log{aiZ^ + (3i)] < oo . (4.11) 

Then the following statements hold. 

(1) The stationary solution (cr^) to (j4.9[) is regularly varying with index a and 
strongly mixing with geometric rate. In particular, there exists a constant 
ci > 0, given in (|4.5|) with Ai ~ aiZ^ + Pi and Bi = ao such that 

P(cr^ > x) - ci X"" . (4.12) 

For Pi = 0, the squared ARCH(l) process {Xf) is regularly varying with 
index a and strongly mixing with geometric rate. In particular, there exists 
a constant cq > given in (|4.5p with Ai = aiZf and Bi — a^Zf such 
that 

P(X2 >x)-cox-". 

(2) Assume a £ (0, 1) U (1, 2) in the GARCH(1, 1) case. Then 

f a e (0,1) 



where h,, 



nE{a^) ae(l,2), 



and the a-stable random variable Z^ has characteristic junction ijja{t) 
exp(-|t|«Xa(i,c+,0)), where 

OO t 

c+ =E[(l+Too)"-T^] and T^ = Y.Il^'^i ^^ + M ■ 

t=i i=i 

(3) Assume a € (0, 1) U (1,2) in the ARCH(l) case. Then 

f a e (0,1) 



(con)-i/" [^X^ - bn) 4 Z„ , where h 



7iE(cr-^) ae(l,2)^ 
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and the a-stable random variable Za has characteristic function ipait) 
exp(-|i|"xa(i,c+,0)), where 



t 

2 



c+ = E[(l + Too)" - T^] and T^ = J2 IT ^'^ 



t=i i=i 



The limit results above require that we know the constants ci and cq ap- 
pearing in the tails of cr^ and X^. For example, in the ARCH(l) case, 

E[(ao + ai ^2)" - (ai X^)"] ElZP" 
Co = - 



a E[(aiZ2)"log(aiZ2)] 



Moreover, (8.66) in Embrechts et al. (|23 ) yields that X"^ = {ao/ai)Too ■ Hence 
the constant cq can be written in the form 

<E[(l + foo)"-r^]E|Z|2" <c+E|Z|2" 

Co - 



a E[(ai Z2)" log(aiZ2)] a E[(ai Z^)°' log(aiZ2)] ' 

The moments of \Z\ can be evaluated by numerical methods given that one 
assumes that Z has a tractable Lebesgue density, such as the standard normal 
or student densities. Using similar numerical techniques, the value a can be 
derived from (|4TT|) . The evaluation of the quantity E[(ao +ai X'^)" - (ai ^2)"] 
is a hard problem; Monte-Carlo simulation of the ARCH(l) process is an option. 
In the general GARCH(1, 1) case, similar remarks apply to the constants c+ and 
ci appearing in the stable limits of the partial sum processes of {cr^). Various 
other financial time series models fit into the framework of stochastic recurrence 
equations, such as the AGARCH and EGARCH models; see e.g. the treatment 



in Straumann and Mikosch (53|) and the lecture notes by Straumann ([55 

In what follows, we consider the GARCII(1, 1) case and prove stable limits 
for the partial sums of (Xt) and (^2). 

Conditions (RV), (MX) and (lACl) 

In what follows, we assume that Z is symmetric, has a positive Lebesgue density 
on M and there exists a > such that (|4.1ip holds. Under these assumptions, it 



follows from Mikosch and Starica (|43|) that (Xf) is regularlyvarying with index 
a and strongly mixing with geometric rate. By Breiman's (|8|) result we have in 
particular, 

P(X2 > x) - E|Z|2" P((t2 > x) - E|Z|2" Ci . 

By definition of multivariate regular variation, the sequence (lATtj) inherits reg- 
ular variation with index 2a from (Xf). By symmetry of Z the sequences 
(sign(Zt)) and (|^t|), hence (sign(Ar4)) and (|A't|), are independent. Then an 
application of the multivariate Breiman result in Basrak et al. (jBl) shows that 
(Xf) is regularly varying with index 2a and 

2a 



V{X >x) = 0.5P(|X| > x) - 0.5E|Z|2" ci X 
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Thus both sequences {Xt) and {Xf) are regularly varying with indices 2a and a, 
respectively. Moreover, (MX) is satisfied for both sequences and we may choose 
fn — rC^ , nin = rC'^ for sufficiently small < 71 < 72 < 1. An application of 
Lemma [XT] yields (|AC[) . Wc omit details. 



4.4.1 Condition ((TB|) for the squared GARCH(1, 1) process 

Recall the notation At = ai Zf_i + Bt = ao, Ilf = n*=i ^^'^ t^^^* i^n) 
satisfies nP{X'^ > a„) ~ 1. The same arguments as for (|4.6p yield 

H h = Zfaf^ \- Zjaj 

= al {Zf ni + --- + Zjlla)+Rd- 

Under the assumption (|4.1ip . E(i?^) < 00, hence P{Rd > = o{P{X^ > a„)). 
This fact and Breiman's result ^ ensure that 

nF{Xf + ---+X^>an) - nP{al{Zfn, + --- + ZjUd)>an) 



E\Z 



2a 



lnnx'>a .^ nz!Tl, + --- + Zjur 

Eizfni + .-. + z^ndi" 



E|Z|2" 

\z^ + zfn^ + ■ ■ ■ + zj_,nd-i\" 



E\Z 



2a 



= 6+(rf). 



In the last step wc used that ^1 is independent of Zi, . . . , and that KA" = 1. 
Write 

Td = Z^.Hi + ■ ■ ■ + Z^^Ud . 

Observe that 

T<iS^arMn2 + --- + n<j+i] . 

The same argument as in Section proves that the right-hand side converges 
a.s. to a finite limit. Hence Td Too for some finite limit T^o = X^t^i n* • 
If a € (0, 1], we have by concavity of the function f{x) = x", x > 0, 

nT^+i - T^] = + Tdr - T^] < < oo . 
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If a e (1, 2) we have by the mean value theorem for some G (0, Zq) and using 
the concavity of the function f{x) = x°'~^, x > 0, 



^ a 



Q-ll 



aE(|Z|2(«-i)) [1 + E(A"^^) + • • • + {E{A"-^)f] 

= aE(|Z|2(«-i))(l-EA°-i)-i < C50. 

An apphcation of Lebesgue dominated convergence yields in the general case 
a e (0, 2) that 

c+ ^ lim[6+(rf+l)-6+(d)]= lim 5KM±^j);^ 



d— )-oo d— >oo E|^| 

^" ' e (o,oo). 



EIZP" 



4.4.2 Condition (|TB|) for the GARCH(1,1) process 

Next we calculate the corresponding value c+ for the GARCH(1,1) sequence 
(Xt). By the assumed symmetry of Z, we have = c_. Slightly abusing 
notation, we use the same symbols b±{d), (a„), Td, etc., as for (Xf). We choose 
(a„) such that nPdA'l > a„) ~ 1. We have 

Sd = Ziai + ■ ■ ■ + ZdCTd ■ 

Since E|Zp" < od we have for any e > 0, 

nY{\\{Zi<Ju...,Zdcyd) - [Zi^i, Z2\A2\°-''<Ju ■ ■ ■ . Zd\Ad\''-''<Jd-i\\ > ea„) 

1/2 



Hence (see Jessen and Mikosch (l36|), Lemma 3.12) 

nF{Sd > a„) - nP(ZiCTi + ZaA^-Vj + ZgAg-'Va • • • + Z^yl^-^ad^i > a„) . 

Proceeding by induction, using the same argument as above and in addition 
Breiman's result, and writing n^^t = Y[l=s s ^ t and list = 1 for s > t, 



28 



we see that 

nV{Sd>a„) - nP(ai(Zi+A0-5Z2 + ---+n°:5Zd)>a„) 

^ [n P(ai > a„)] ]E[(Zi + A^'^ Z2 + • • • + U°-}Zd)f] 

E[(Zi + Z2 + • • • + nl^^Zd)l"] 



E|Z|2« 

E[|Zi+A0-5Z2 + ---+n0-5z<i|2«] 



= fe+(d). 



2E|Z|2" 

In the last step we used the symmetry of the Zt's. Writing Tj, ~ X]f=i ^i^2'i' j 
we have 

and 

00 00 

i=l i=l 

Since E[(logA)i/2] = 0.5ElogA < 0, the right-hand side converges a.s. to a 
finite limit. By a Cauchy sequence argument, 

00 

for some a.s. finite Too. 

Let be an independent copy of (^1,^2), independent of Td. Assume 

2a G (0,1]. Then by symmetry of the Zt's and since E[(yl2)"] — 1, using the 
concavity of the function f{x) = x'^", a; > 0, 

E[||Td+i|2"_|Tdp"|] 

= E[\\Z[ + (A^)0-5T,|2" - 5t,P"|] 

= E[\\Z[\ + {AX' Tdl^'' ~ liAX'Tdn] 

= E[{\Z[\ + {AX' {Td)+f'' - {{AX-HTci)+f"] + 

m{\Z{\ - {AX-' {Td)-)"' ~ ((A^)"-5(T,)_)2"|/{(^,)o.= (^,)_^|z;|}] + 
mHAX-' (Td)- l^il)'" - {{AX-HT,)-)"'\Iua',)o.ht,)_>iz',\}] 

s; E|zp". 

For 2a G (1,2) we use the same decomposition as above and the mean value 
theorem to obtain 

E[||Td+i|2"-|rd|2"|] <; 2aE[{AX-'\Td\ + \Z{\]^'''^ 

^ 2aE[{AX-'\Td\?°'~^ + 2aE|Zp"-i . 
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The right-hand side is bounded since EjZp" < oo and, using (pi^ and E[^"-°-5] < 
1, 

OO OO 

1=1 i=l 

Now we may apply Lebesguc dominated convergence to conclude that the limit 



c+ = lim [b+{d+l) - b+{d)] 



2E|Z|2" 



exists and is finite. 

Since Too and Z( assume positive and negative values we have to show that 
c+ > 0. First we observe that 



c+ = lim rf-iE[|rrf|2«] , 



(4.14) 



Applying Khintchine's inequality (see Ledoux and Talagrand (|4l[ )) conditionally 
on (1^(1), we obtain for some constant > 0, all 1, 



n\Tc 



|2qi 



1=1 

d 

= c„«-"E(^(aiZ.f +/3i)n2.,-/3i;^n2 

4=1 1=1 

r"[E(5]n2,.+i)"-/3rE(^n. 

i=l 

r"(i-/3nE(En2.) 



i=l 



Now, Remark lOl and (|4.14p imply that c+ > 0. 



4.4.3 Condition 5 

We assume Z symmetric. Then EX = for 2a G (1,2) and (jCTI) holds for (Xt). 
For {X'^), (jCTp cannot be satisfied and needs special treatment. If a G (1,2) 
all arguments above remain valid when Xf is replaced by E{X^) — E{X'^). 
We summarize our results for the GARCH(1, 1) process (Xt) and its squares. 

Proposition 4.8. Let (Xt) be a strictly stationary GARCH(1, 1) process with 
symmetric iid unit variance noise (Zt). Assume that Z has a positive density 
on R and that (|4.1ip holds for some positive a. Then the following statements 
hold. 
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(1) The sequences {Xt) and (Xf) are regularly varying with indices 2a and a, 
respectively, and both are strongly mixing with geometric rate. In particu- 
lar, 

F{X > x) - ^E|Zp" cia;-2" , 

where c\ is defined in (I4.12p . 

(2) ^sswme 2q; e (0,2). Then 

where Z2a is symmetric 2a-stable with characteristic function 4'2a{t) = 
exp(-|tp"X2a(i,c+,c+)), where 

^ E[|Zo + laiZg + /3i|°-5roop" - llaiZg + /3i|"-5Toop"] 

2E|Z|2" 

and 

00 t—i 



t=l i=l 

(3) ^ssitme a G (0, 1) U (1, 2). Then 

n 

(Ci i?|Z|2°n)-l/" {Y.^t- bn) A ^a, 
t=l 

where 

Jo if a e (0,1) 

" ~ \ nE(X2) if ae (1,2), 

and is a-stable with characteristic function 

V'c.(i) =exp(-|irxa(i,c+,0)), 

where 

^ E[(zg + f^)"~r^] 

+ E|Z|2" 

/n t/ie above relation, 



t=l i=l 



Remark 4.9. For ARCH(l) processes the above technique of identification of 
parameters of the hmiting law was developed in (Q). 
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4.5 Stable stationary sequence 



In this section we consider a strictly stationary symmetric a-stable (sas) se- 
quence {Xt), a € (0,2), having the integral representation 

Xn^ f fnix)M{dx), neZ. 

JE 

Here M is an sas random measure with control measure /i on the cr-field £ on E 
and (/„) is a suitable sequence of deterministic functions /„ S L" {E , £ , fj.) . We 
refer to Samorodnitsky and Taqqu ((sol ) for an encyclopedic treatment of stable 
processes and to Rosinski (|49| ) for characterizing the classes of stationary (Xt) 
in terms of their integral representations. 
Then for some sas random variable Ya, 



E 



{h{x) + ---+Ux))M{dx)^Y^([ \h{x) + ■■■ + UxT^i{dx)V'" . 

J E 

(4.15) 

Since P{Ya > x) ^ Q.^cqx " for some cq > (see Feller (|26l )). we have with 
nP{\X\ > an) ~ 1, 

""^^^'^""^ /,|/i(x)hMrf^) 

Moreover, it follows from (|4l^ that a-^Sn 4 for some Za if and only if 

n-^b+{n) c+ (4.16) 

for some constant c+ and the limit Z^ is sas, possibly zero. 

Since we know the distribution of a~^Sn for every fixed n we do not need 
Theorem 13.11 to determine a sas limit. In the examples considered above we 
are not in this fortunate situation. In the sas case we will investigate which 
of the conditions in Theorem 13.11 are satisfied in order to see how restrictive 
they are. Since the finite-dimensional distributions of {Xt) are a-stable, (RV) 
is satisfied. Conditions (|CT[) and EX = for a G (1,2) are automatic. Under 
(|4.16p . using the special form of the characteristic function of a sas random 
variable, (MX) holds for any sequence m„ oo. Condition (|AC[) is difficult 
to be checked. In particular, it does not seem to be known when (Xt) is strongly 
mixing. An inspection of the proof of Lemma 13.41 using the particular form of 
the characteristic functions of the sas random variables, shows that ()ACp can 
be replaced by (|TBp which implies (|4.16p . Thus (|TBp is the only additional 
restriction in this case. 
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